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Abstract 

We analytically investigate the kinetic Gaussian model and the one- 
dimensional kinetic Ising model on two typical small- world networks (SWN), 
the adding-type and the rewiring-type. The general approaches and some ba- 
sic equations are systematically formulated. The rigorous investigation of the 
Glauber-type kinetic Gaussian model shows the mean-field-like global influ- 
ence on the dynamic evolution of the individual spins. Accordingly a simpli- 
fied method is presented and tested, and believed to be a good choice for the 
mean-field transition widely (in fact, without exception so far) observed on 
SWN. It yields the evolving equation of the Kawasaki-type Gaussian model. 
In the one-dimensional Ising model, the p-dependence of the critical point is 
analytically obtained and the inexistence of such a threshold p c , for a finite 
temperature transition, is confirmed. The static critical exponents, 7 and (5 
are in accordance with the results of the recent Monte Carlo simulations, and 
also with the mean-field critical behavior of the system. We also prove that 
the SWN effect does not change the dynamic critical exponent, z = 2, for this 
model. The observed influence of the long-range randomness on the critical 
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point indicates two obviously different hidden mechanisms. 
PACS number(s): 89.75.-k, 64.60.Ht, 64.60.Cn, 64.60.Fr 
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I. INTRODUCTION 



Small world networks (SWN) are those intermediate between a regular lattice and a 
random graph (see and references therein for review). They are believed to catch 

the essence of networks in reality, such as neural networks 0], power grids, social networks, 
and documents on World Wide Web |],[7[], where remote vertices, while locally clustered, 
often have the chance to be connected via shortcuts. Since 1998, when Watts and Strogatz 
presented a simple model showing SWN effects it has been studied intensively and 
extensively M. From the point of view of statistical physics, the presence of shortcuts 
assists the system to behave as a whole, showing global coherence and new, SWN behavior, 
possibly apart from their ordinary properties people are familiar with. 

One may be curious about, to what extent the features of phase transitions will be dif- 
ferent in spin-lattice models built on small-world networks. Because the SWN effect widely 
exists in reality, this question is also of much significance. Although new and interesting 



features have been revealed recently pL|, 10, 11, 12, 13, 14, 15, 16], it is still far from being com- 



pletely answered. The dynamic aspect has been even less well understood ||17|| . Naturally 
we may expect the evolution of a single spin to be influenced partly by the overall system, 
however we find it difficult to offer more specific information. This is due to the complexity 
of the dynamics itself, and the often formidable mathematical task. 

In this article, we report our work in this interesting problem — the critical dynamics 
of spin-lattice models built on SWN. As the introductory content we shall first discuss the 
general approach. Then, among the various model systems, we choose two special ones 
for a detailed investigation: the Gaussian model, which is relatively easy in mathematics, 
and the one-dimensional Ising model. New and interesting kinetic features are revealed 
analytically and the study in the dynamic aspect also yields much information about the 
static properties, such as the critical point. 

This article is organized as follows: In Sec. |TJ we give the definition of spin-lattice models 
built on SWN. Sec. [TTI| contains the discussion of the general approach (along with a brief 
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review of the dynamic mechanisms). The direct application on the kinetic Gaussian model 
can be found in Sec. [TV} In Sec. |V|, we present a simplified method and prove its validity in 
Sec. |V A| by comparing its result with the rigorous one obtained in |V[ Sec. |V B| and Sec. 
V C| are devoted to the further applications of this method on the Kawasaki-type Gaussian 
model and the one- dimensional Glauber-type Ising model, respectively. In Sec. [VT], the 
influence of the randomness on the critical point is analyzed. Sec. |VII| is the summarization 
with some discussions. 



II. SPIN-LATTICE MODELS BUILT ON SMALL- WORLD NETWORKS 

Following the first prototype IJ of SWN, there have been a number of variants (see, 
for example, Ref. |18|]) in two basic groups, which can be constructed as the following: the 
initial network is, for example, a one dimensional loop of iV vertices, each vertex being 
connected to its 2k nearest neighbors. (1) Each pair of random vertices are additionally 
connected with probability pa\ (2) the vertices are then visited one after the other, and each 
link connecting a vertex to one of its k nearest neighbors in the clock wise sense is left in 
place with probability 1 — p^, and with probability p^ is reconnected to a randomly chosen 
other vertex. We may call the first group adding-type small-world networks (A-SWN) and 
the second group rewiring-type (R-SWN). Both of these modifications introduce long range 
connections. The above algorithms can be extended to systems with higher dimensionality 
and even fractal structures. (For example, for a spin located somewhere in a cubic lattice, 3 
of the 6 bonds in some fixed directions will be reconnected with probability p^ respectively.) 
These constructions, with any (infinitesimally low) fraction of shortcuts, allow to reconcile 
local properties of a regular network (clustering effect) with global properties of a random 
one (the average distance I ~ log 10 N) . 

A spin-lattice model built on SWN, which is the focus of this article, can be thus de- 
fined: In a .D-dimensional regular network consisting of N spins with periodic boundary 
condition, each spin is linked to its 2kD nearest neighbors (in this article we will choose 
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k = 1). Then, (1) a certain number of supplemental links are added; (2) a portion of the 
bonds are rewired. Whether it is in group (1) (A-SWN) or group (2) (R-SWN), two spins 
connected by a shortcut act in the same way as those connected by a regular bond, that 
is, their interaction contributes to the system Hamiltonian and they obey the redistribution 
(exchange) mechanism. In short, there is no difference between a long range bond and a 
regular one. Other definitions may exist. For example, because these small-world networks 
are model systems for the networks in reality, one may have reason to choose certain phys- 
ical quantities, such as J in the Ising model, to be different. However, although it is not 
necessary, we will limit the scope of this article to the above rule. 



III. THE DYNAMIC MECHANISM 

The various dynamic processes in the critical phenomena are believed to be governed by 
two basic mechanisms, Glauber-type with order parameter nonconserved and Kawasaki-type 
with order parameter conserved. Recently, a Glauber-type single-spin transition mechanism 
T9| , |20|| and a Kawasaki-type spin-pair redistribution mechanism have been presented 



as the natural generalizations of Glauber's flipping mechanism [^] and Kawasaki's exchange 
mechanism [2J], respectively. They are generally applicable and mathematically well orga- 



nized. Our work begins with a review of these two mechanisms. It is for your better 
understanding of the calculations in Sec. |IV| and [V], and it also might be a convenient 



reference for further studies. 



A. Single-spin transition mechanism 

Glauber's single-spin flipping mechanism allows an Ising system to evolve with its spins 



flipping to their opposite. In the single-spin transition mechanism [|1^], a single spin a, may 
change itself to any possible values, <7j, and the master equation is 

j f P(W}, *) = -££ [Wifo - &i)P({*}, t) - Wtfi - a h t)] . (1) 



The transition probability is in a normalized form determined by a heat Boltzmann factor, 

1 



Wi(ai -> <Ji) = — exp [-(3H ({crjjn} , , <3i = £ ex P [~/^ > o-j)] 

where /3 = 1/ksT. Based on the master equation, Eq. ([!]), one can prove that, 
dq k (t) 



(2) 



dt 



-Ik (t) + £ 
W 



P(W;t) 



(3) 



where (£) = S{o-} ojfef ({°"} j There are also corresponding equations for the correlation 
functions. 

B. Spin-pair redistribution mechanism 



Kawasaki's spin-pair exchange mechanism allows an Ising system to evolve with its near- 



est neighbors exchanging their spin values. In the spin-pair redistribution mechanism [21 



two connected spins, <jj and 07, may change to any possible values, 5j and 07, as long as 
their sum are conserved. The master equation is 
d 



dt 



P(W}, *) = ££ [-W fl (* j( Ti - P (M ; t) 

(jl) 

+ Wjl (frjfrl -> OjOl) P ({(Tift, CT/^fc} , CTj, erf, t)] . 



(4) 



The redistribution probability is also in a normalized form determined by a heat Boltzmann 
factor, 



Wji (<Tj<ri -> = ^-^+0,,^+^ exp ({cr m } m ^ y , o-j, a; 

where the normalization factor Qji is 

Qji = £ s ^+a l ,a J +a l exp {{^m} m ^j,i > ^ 



(5) 



Based on the master equation, Eq. ([§), one can prove that, 
dqk (t) 



dt 



-2Dq k (*) + ££ 

{a} W 



£ &kWk,k+w {.&k&k 



P(W;t), (6) 



where .D is the dimensionality and stands for the summation taken over the nearest 
neighbors. 



C. How to apply them on SWN 



In the construction of SWN, according to a certain probability, we will have a whole 
set of possible realizations. So the theoretically correct way of treating dynamic systems 
built on SWN actually consists of three steps: First we have to make a full list of all the 
possible realizations and point out the probability of each one of them. Second, we treat 
each system respectively (apply the dynamic mechanism and obtain the master equation 
and the physical quantities of interest). Third, we give the expectation value with all these 
results. This is cumbersome, but conceptually straightforward. In Sec. [V| we shall discuss 
whether there is any simplified method. (In fact, there is.) 

IV. KINETIC GAUSSIAN MODEL GOVERNED BY THE GLAUBER-TYPE 

MECHANISM 



We have discussed the general approach right above in Sec. [Ill Q, and in this section we 



will directly perform the calculations according to the three steps, to study the 3D kinetic 
Gaussian model governed by the Glauber-type single-spin transition mechanism (see Sec. 
pi A|) . The calculation to be carried out is very long but fortunately we can borrow some 
results from our earlier studies [plpOf . 

The Gaussian model, proposed by Berlin and Kac, is a continuous-spin model. Its 
Hamiltonian, 

-/3H = Kj2°i°i, (7) 

(id) 

where K = J/ksT. The spins can take any real value between (—00, +00), and the 
probability of finding a given spin between and a^+dak is assumed to be the Gaussian- type 
distribution, / (cr*.) do^ ~ exp (— \ da^, where b is a distribution constant independent 
of temperature. Thus, the summation for the spin value turns into the integration ~ * 
f^oc f ( (J ) dcr. This model has been studied often as a starting point for investigations of 
other systems. 



Governed by the Glauber-type single-spin transition mechanism, the expectation value 
of single spin obeys [|19[ 

dqk(t) /+\, K u\ fa\ 

—jj— = ~q k W + y L Qv (*) • ( 8 ) 

(<7 fc '><?fc) 

Since on SWN a spin may have a neighbor located very far from it, in the following calcu- 
lations the summation is to include every spin that is connected with a^. Taking average of 
Eq. (H), we can obtain the evolution of the magnetization M (t) =J2kQk (t) /N. 

A. kinetic Gaussian model on an adding-type small-world network 

First we treat the 3D kinetic Glauber-type Gaussian model on A-SWN consisting of N 
(a very large number) spins. The Kawasaki-type is still tractable, but the mathematical task 
will be more complex. We will leave that till Sec. 0. Besides the regular bonds, each pair 
of spins, no matter how far apart, is connected via an additional bond with probability p^. 
Actually there are 2 N ( N ~ 1 ^ 2 different networks, each with a given probability. As a model 
system for the networks in reality, we expect the number of these bonds, n ~ N (N — 1) pa/2, 
to be much smaller than N. So practically we require PaN <C 1. With respect to a specific 
spin aijk, all the networks can be divided into N groups listed as the following. 

(0). There is no random bond on cr^, and the probability is (1 — Pa) N ~ 1 - According to 
Eq.® 

d K 

^lijk if) = —Qijk W + yX! (<li+w,j,k (t) + Qi,j+w,k (t) + qi,j,k+w (t)) ■ 



(1). There is only one random bond on cr^, and the probability is C^^pa (1 _ Pa 



N-2 



In fact this group can be further divided into N — 1 subgroups, each corresponding to a spe- 
cific spin connected to cr^, and each with the same probability p^(l — Pa) N ~ 2 ■ Averaging 
them we get 

~|~ = ~9ijk (*) + y (li+w,j,k (t) + qi,j+w,k (t) + qi,j,k+w (*)) + (t) . 



s 



(n). There are n random bonds on a^k, and the probability is C]^_ 1 p^ (1 — Pa) 
Similarly, in average we get, 

= -lijk (*) + y £ (ft+wj,* (*) + ftj+w,* (*) + ft,j,fc+t» (*)) + ~y~ M (*) • 



N-n-1 



(N — 1). There are iV — 1 random bonds on er^, and the probability is p A 

dqij ^ = -%k W + yE (9i+wj,k (t) + Qi,j+ W ,k (t) + %j,k+ w (t)) + j(N-l)M (t) 
Thus, over all the realizations, 

" = ~Qijk (t) + — ^2 {qi+ w ,j,k (t) + qi,j+ w ,k (t) + qi,j,k+w (t)) 



K 



N-l 



+ 1 -HC n N . 1 p n A {l-PAf- n ~ 1 nM. 



(9) 



n=0 



Fortunately we find the following relationship 



N-l 



E Cft-iJfl (1 - PA)"-"- 1 n=(N-l) p A , 

n=0 



(10) 



and thus 

= (*) + y £ (ft+wJ.* (*) + Qij+w,k (*) + ?ij,k-Ho (*)) + y ~ !) PaM (t) , 

(11) 

Similar results can be found in one and two dimensional models, and taking average we 
obtain 



dM (t) 
dt 



2DK K . 

i— i— j-(w-i)p* 



M{t) 



(12) 



where D is the dimensionality. The solution of Eq.(|T2"D is 



M (t) = M (0) exp (-t/r) 



where the relaxation time 



and the critical point 

K * = 2D + (N-l) PA (13) 
It is well known that the critical point of the regular model K r c e9 = b/2D, we can clearly see 
that the critical temperature will get higher as more long range bonds are added. Actually, 
for a vertex located on an A-SWN, the number of the long range bonds ~ (N — 1)pa- 
In the small-world region, where the expected ha for most of the vertices is very small, the 
change of the critical point will be almost unperceivable. (The analysis of this result can be 



found in Sec. VI. 



B. kinetic Gaussian model on a rewiring-type small-world network 

Second we treat the kinetic Gaussian model on a R-SWN with characteristic probability 
Pr. Because of the length, here we only give the details of the ID case. With regards to a 
specific spin a k , all the networks can be divided into four major groups listed as the following, 
(1) both of the two regular bonds on a k , connecting o~k-i and (Xfc+i, are not rewired; (2) the 
bond connecting a^-i and a k is left unchanged but that connecting a k and ak+i is rewired; 
(3) the bond connecting o~ k and a k +i is left unchanged but that connecting c^-i and a k is 
rewired; (4) both of them are rewired. The final result comes from the summation of all the 
four parts (see Appendix A for the details): 

^1 = _ qk {t) + p J^M (t) + (1 - p R ) j (q k ^ (t) + q k+l (t)) , (14) 

In two and three dimensional models, there are 2 4 and 2 6 major groups respectively, and 
they can be treated in the same way. Similar results are obtained: 



dM (t) / 2DK\ 



dt ' ' 

Obviously the relaxation time 
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-)M(t). (15) 



T 1-K/K*' 

where the critical point 

K? = K T c eg = b/2D. (16) 

To summarize, we strictly obtain the evolution of the Gaussian model built on SWN. On 
A-SWN, the critical temperature will get higher as more long range bonds are added, while 
on R-SWN the critical temperature is unchanged. With the dynamic scaling hypothesis 
t ~ <^ z ~ \T — T c \~ zv , we have zv = 1, where z is the dynamic critical exponent and v is 
the correlation length critical exponent. We shall leave the summarization of its properties 
to Sec. |V11[ The influence of the long range bonds on the critical point, which is a rather 



interesting topic, will be discussed in Sec. |VT| . 

The Gaussian model, being an idealization, often has the value of serving as a starting 
point for the more general studies. The evolving equations of the individual spins show 
distinctly the influence of the global coherence, which automatically takes a mean-field-like 
form as if it comes from an averaged spin. This character, as well as some already-proved 
facts in earlier studies about the static behavior, provides us with some helpful hints for the 
simplification of the method, which is the topic of the following section. 

V. THE SIMPLIFIED METHOD 

Presently, for the critical dynamics on SWN, we still lack a well-established approach, 
which should be both theoretically reliable and practically feasible. The mentioned-above 
method (and the results) is theoretically rigorous, but it is too complex for the other models, 
even the simplest one-dimensional Ising model. In the above section, we have pointed out 
that, in the dynamic evolution of the Gaussian model built on SWN, the influence of the 
system as a whole on individual spins is mean-field like. On the other hand, recent studies on 
the Ising model and the XY model built on SWN have also shown that the phase transition 
is of the mean-field-type (see 0|^,[l^,[ll],[ll|,[l5|] and Sec. [V Q for details). The various model 
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systems built on SWN probably belong to the same mean-field universality class, and they 
might be treated with the same dynamic approach. In this problem, which is found to be 
of the mean-field nature, the following method may be the best choice. We deem all the 
possible networks as a single one. The effective Hamiltonian of a spin-lattice model built on 
such a network is defined as the expectation value over all possible realizations. Its effective 
behavior, e.g. the redistribution between two vertices connected with each other, is also the 
averaged result. For example, two sites, Oi and &j, are connected with probability p, then 
redistribution occurs between them with probability pWij. (It is the basic assumption when 
we are treating a system governed by the Kawasaki- type mechanism) . Then we can directly 
apply the dynamic mechanisms to this system. This method, which is of the mean-field 
nature, is certainly more tractable in mathematics. We apply it to the kinetic Gaussian 
model again and see if it will lead to the same result. 



A. Application in kinetic Gaussian model governed by the Glauber-type mechanism 

In the 3D Gaussian model (N spins in total) built on A-SWN with periodic boundary 
condition, the effective Hamiltonian 



— (37i — K (Jijk (<Ti+i,j,k + &i,j+i,k + 

i,j,k 

I i 

+ 2 KpA ^ ° { i k ^ ai 'i' k ' ~ 2 KPA ^ °* jk ' 

~ ij',fc i',j',k' ~ ij,k 



(17) 



With the Glauber-type transition mechanism, we begin to derive the single-spin evolving 
equation according to Eq.(H), along the same line of the calculations in Ref. |nj. We obtain 



K 

^ ^ijkWijk {&ijk —* &ijk) — ~g 



^ (Vi+w,jk + &ij+w,k + &ij,k+w) + (N — 1)pa& 



(18) 



where 



a 



\lmn / Imn 

Substituting Eq. (|18D into Eq. fl|) we S et 



&lr 
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<]<I " U - -q ijk (t) + T E (?*+«jfc (t) + Qij+ W ,k (t) + q ih k+ W (*)) + x - !) PaM (t) 



at U)=±1 



(19) 



Similar results can be obtained for one and two dimensional models. These results, along 
with those for M (t), are exactly in accordance with what we have obtained in Sec. B 
Eqs. (0) and ©. 



In the 3D kinetic Gaussian model built on R-SWN, the effective Hamiltonian, 

— (3H — K (1 — p R ) Y a ijk (&i+l,jk + &i,j+l,k + ^ij.fc+l) 

3 3 
+ K PRJjJ2 a ijk a i'3'k< - Kp R —Y. a lk- (20) 

i,j,k i',j',k' i,j,k 

Similar calculations yield 

dq^k 6 A" K y-^ . 

^ = — ^ijfc + PR—^-M + (1 — PR) ~g {qi+w,j,k + qi,j+w,k + qi,j,k+wj ■ 



This is also in accordance with the rigorous result, Eqs. (|14] , |15| ) . Thus, in the kinetic 
Gaussian model this simplified method yields the same results as the rigorous ones obtained 
with the more complex standard method. As various spin-lattice models, such as the Ising 
model, are believed to show mean-field behavior on S WN, we believe that this simplified 
method is able to provide at least qualitatively correct information. In this sense, it is 
very different from the mean-field approximations taken in other universality classes. In 
later studies, we will use it to study some more complex problems. The following are two 
examples. 



B. Application in kinetic Gaussian model governed by a Kawasaki-type mechanism 

Now we apply this simplified method to study the diffusion process in the kinetic Gaus- 
sian model built on S WN. Although we still have to deal with many complex equations, it 
is relatively easy compared with the formidable task of the standard approach. In such pro- 
cesses, the system is governed by the Kawasaki-type redistribution mechanism. As already 
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mentioned at the beginning of this section, the system behavior, just as the Hamiltonian, is 
also averaged over all possible realizations. The basic equations for the Kawasaki-type dy- 
namics listed below are generally applicable in various order-parameter-conserved processes. 
(1) On .D- dimensional A-SWN: Accordingly the master equation should be modified as, 
d 



dt 



+ Wjl (d-j&l -> CTjUl) P ({(TijLj, ai^k} , <7j, en; t)] 

+\pa E E E [-Wji fan - *&) P (M ; t) 
+Wji fa&i — > OjOi) P {{a^j, <Jij=k) , (tj, o-/; t)] , 



(21) 



where the redistribution probability Wji fan — y &jn) is of the same form as Eq. @. With 
Eq. (pi]) we can get that 
dq k (t) 



dt 



-2Dq k (t) + E E 

M w 



E °kW k;k 

+w 



+wj 



+Pa \-(N -l)q k (t) + E 



P(M;t) 
P(W;t) 



(22) 



(2) On D-dimensional R-SWN: Accordingly the master equation should be modified as, 
d 



dt 



P({a}, t) = (1 - p fl ) E E fan - &&) P (fa} ; t) 

+Wji faa t — >■ ajai) P ({<JijLj, n^k} , cr j7 at; t)] 
+DpR N^l E E E [-Wji fan - P (fa} ; t) 
+Wji faai — >■ ajai) P ({n^j, n^k} , a j: a t ; t)] . 



The redistribution probability Wji fan — > <7j<5"z) is of the same form as Eq.(|5], and 



(23) 



dt 



{l-p R ){-2Dq k (t) + EE 



E 

,k+w i&k&k+w > &k&k+'b 



P({a};t) 



Dp 



N 



\\-(N-l)qk(t)+Y: 
{ M 



E E °kW k i fan — > a k ai) 

l^k a k ,ai 



P(fa};t) 



[i-pr)R 



(1) . ^Pil „(2) 



N — 1 



(24) 
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(1 2) (1 2) 

Although A\ ' ' and R k ' ' are of the same form respectively, we use different symbols be- 
cause they are actually different (determined by the Hamiltonian-dependent redistribution 
probability W). 

First we treat the Gaussian model built on A-SWN. The first part, A\j k , comes from the 
redistribution between the nearest neighbors. In the 3D case, with Eq. ([|) we can obtain 

A ifk = 2(b + K) b ^ qi+1 ' j ' k ~ qijk ^ ~ ( qijk ~ qi ~ 1J ' k ^ 

+ — lijk) — (Qijk — Qi,j-l,k)} + [{Qi,j,k+1 ~ Qijk) ~ (Qijk ~ Qi,j,k-l)]} 

K 

+ 7777 — ; — t77 (29i-l,j,fc — Qi-l,j+l,k ~ Qi-l,j-l,k) + {^Qi-l,j,k ~ Qijk ~ Qi-2,j,k) 
I [p + K ) 

+2 (2gj + ij ) jt — gj+ij+i^ — gi+i.j-^fc) + (2qi + xj i k — q^k — qi+2,j,k) 
+2 (2gj J _ 1) fc — gi^-— ljfe+i — qij-i^-i) + (2<li,j~i,k — Qijk — Qij-2,k) 
+2 (2gj J - + i i fc — qij + i 7 k + x — qij+i^-i) + (2(^+1^ — q^ — qij+2,k) 
+2 (2qij t k-i — Qi+i,j,k~i ~ Qi-ij,k-i) + (2gij,fc-i — — qij±-2) 

+2 (2qij i k+l — Qi+l,j,k+l — Qi-l,j,k+l) + (2gij,fc-fl — — gjj,A:+2)] • (25) 

Actually this is not as complex as it seems. With the lattice constant a we can transform 
the above expression to be 

A % = rn<(l- K ) v2 '>^ t) - (26) 

(2) 

The second part, A\-' k) comes from the redistribution between and the farther spins, 
<Ji>j>k>- The result is 

Af k = -(N -\)q t]k (t)+ l - (qijk(t) + qi>fk>(t)) 

i'j'k'y^ijk 

K 

+ 77T ^ ^ [(qi+w,j,k (0 + Qi,j+w,k (0 + qij± +w 

(*)) 

20 i'j'k'^ijk w 

~ (Qi'+w,j',k' CO + Qi',j'+w,k' (0 + Qi',j',k'+w (0)] 5 (27) 



Substituting Eqs.(|^) and (|27|) into Eq . (|22|) , we obtain the evolving equation of the Gaussian 
model built on 3D A-SWN, 



15 



dq (r, t) 
dt 



3a 2 



'b 
6 



b + K 
PA (N-l) 



k)+ Pa {N-r 



a 2 K 
~2b~ 



V 2 g(r,t) 



Q —)(M(t)-q(r,t)). 



(28) 



Similarly, on 3D R-SWN we obtain 



dq (r, t) 
dt 



(1-Pr) 
3pn 



3a 2 



1 



b + K(l-p R ) 
6K 



6 



-K{l-p R ) 



K I 

+ 3p^(l-Pfl) V 2 g(r,t) 



J- -PR) 



(M(t)-q(r,t)) 



(29) 



From these two equations we can clearly see the influence of the system built on SWN 
as a whole on individual spins. On regular lattices the evolution of the system can be 
explained by the diffusion mechanism, while on SWN, to some degree, the individual spins 
will automatically adjust itself to approach the average magnetization. This is obviously 
the result of the global coherence introduced by the small fraction of the long range bonds. 

On regular lattice, where the evolution is pure diffusion, dq (r, t) /dt = W 2 q (r,i), and 
the diffusion coefficient T> will vanish near the critical point. However, on A-SWN, where 
dq (r, t) /dt = V'V 2 q (r,t)+C [M (t) - q (r,t)], two temperatures, T£ and T£ can be obtained 
by setting V and C to be zero, respectively. T c ^ will be lower than the critical temperature on 
a regular lattice, T£ eg , and suggests that the point at which the diffusion stops will be lowered 
by the randomness. T c ^ equals T c res , and at this point the evolution will be pure diffusion. 
Similarly there are also two temperatures for R-SWN, but we will have < < T^ eg . 

Obviously, the system behavior strongly depends on the temperature. For example, we 
study a one- dimensional system and the initial magnetization is q (x, 0) = sinx. When the 
temperature T > T c2 , both T>' and C are positive, and the magnetization will approach 
homogeneity. When T < T cl , both V and C are negative, and the inhomogeneity will 
be getting more remarkable during the evolution. When T c2 < T < T cl , V > but 
C < 0, and this is a more complex region. Although here we can still easily predict the 
system behavior with Eq. (p8~D and obtain a stationary point, generally the evolution will 
be strongly dependent on the local magnetization. 
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C. Application in Ising model 



The second application of the simplified method is on the one-dimensional ferromagnetic 
Ising model governed by the Glauber-type mechanism. Recently it has been studied both 
analytically [fflBI and with Monte Carlo simulations |T0| , p!T| , p!2| , p!3| , p^| . Due to the mathe- 
matical difficulties, the study was not going very smoothly at the beginning. In Ref. ||, 
Gitterman concluded that the random long-range interactions, the number of which above 
a minimal value, lead to a phase transition. In Ref. M, Barrat and Weigt used some ap- 
proximations and expected a finite critical point " at least for sufficiently large p and k > 2" 
on R-SWN. They found out that this transition is of the mean-field type. Although they 
could not calculate the transition analytically, the numerical computation demonstrated a 
nonvanishing order parameter in the presence of a vanishingly small fraction of shortcuts, 
for k = 2 and k = 3. The numerical results seemed to support the following relationship, 
T c ~ —2k/ log 10 pR. The more recent Monte Carlo studies have proved the above-mentioned 
conclusions on A-SWN and R-SWN respectively, but in those cases k — 1. The analysis of 
the relationship between T c and pr can be found in Ref. ||11|| . The critical exponents obtained 
from the simulation JTT|JT2]JT5[| , (3 ~ 1/2, a ~ 0, and v m 1/2, further establish the mean-field 
character. However, although there is substantial numerical proof, the inexistence of such a 
threshold p c , for a finite temperature transition, is still to be confirmed. 

In this article, we will apply the simplified method, which is of the mean-field nature, to 
study the dynamic properties of the one-dimensional Ising model built on both A-SWN and 
R-SWN. We hope the success in the Gaussian model will continue in Ising system, of which 
the behavior is already known to be the mean-field-type. Although we can not obtain the 
full picture of the evolution, we are able to get the exact p-dependence of the critical point, 
and some interesting critical exponents. As will be shown below, our result agrees perfectly 
with the above-mentioned numerical simulation. 

We find that the system shows very similar behavior on A-SWN and R-SWN. We shall 
give the details of R-SWN only, and report the results of A-SWN later. 
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On a R-SWN, the effective Hamiltonian is the same as Eq. fl20|) (the one-dimensional 
version). We substitute it into the single-spin evolving equation, Eq.(||) and obtain, 

= ~q k (t) + E tanh \ K (! - PR) iPk-i + fffc+i) + 2Kp R a] P ({a} ; t) , (30) 

M 

where a = J2j^k °j/ {N - 1) ~ Ej Cj/iV. 

If p# = 0, then it is the one-dimensional Ising model on regular lattice, which we are 
familiar with. One can continue to write 

dq k it) , . 1 , , 

— t; — = -9* (*) + - (9fc-i + ft+i) tanh2K, 
at 2 

and 

= -M(t)(l-tanh2K). 

Civ 

It yields M (t) oc e"* /r , where r = (1 - tanh2K)~ 1 . When K -> if c re9 - oo, r ~ £ z ~ 
(e 2i ^) , and thus z — 2. 

If the rewiring probability p# = 1, then 

= -ft (t) + E ^nh P ({a} ; t) , 

M 

and 

^M^l = - M it) + E ^nh (2Ka) P {{a} ; t) . 
dt {*} 

Although (tanh (2Ka)) ^ tanh (2KM), if we study the case when, near the critical point, the 
system is in almost thorough disorder, we will have (tanh (2Ka)) ~ tanh (2KM) ~ 2KM, 
and 

~ -M (t) + tanh [2KM (*)] . 

CI I 

This helps to determine the critical point, K c = 1/2. When K < K c = 1/2, then the 
system will be stable in a disordered state with M = 0, but when K > K c = 1/2, there 
appears some kind of order. Taking Taylor expansion, one can find that when K is near 



1 /2 

K c , M ~ (K — K c ) 1 . This leads to (3 — 1/2 (in this specific situation). We also get the 
relaxation time r ~ \K — K c \~ x , with the scaling hypothesis (see below) we have zv — 1. 

Qualitatively similar situation should also be found when p R is between and 1. First 
we assume that there is a critical temperature, above which the system is disordered and 
below which there begins to show nonzero magnetization. From Eq.([3(]) we find that if 
initially M = is given, then the system will stay in this disordered state. But below the 
critical temperature this equilibrium will not be stable. We can determine the critical point 
introducing a small perturbation. When M — * 0, 

' : " ~ -Qk (t) + E tanh pr) K-i + vk+i)] p (M ; t) 

M 

+2Kp R 5>{l" tanh 2 [K (1 - p R ) {a k ^ + a k+1 )]} P ({a} ; t) 
M 

= ~Qk (*) + ^ (9fc-i + tann t 2 ^ C 1 _ P«)] 

+2#p B £ a |l - 1 (1 + a k ^a k+1 ) tanh 2 [2tf (1 - p R ))\ P ({a} ; t) . (31) 
M 

Because M is very small, we believe 2 fe (^Ofc-iOfc+i) is an even smaller quantity of higher 
order. Thus 

dM(t) 



dt 



dt 



-M + Mtanh [2K (1 - p fl )] 
+2K^ |l - ~ tanh 2 [2tf (1 - p R )]} M. 



The critical point can be determined as 
tanh 



2K? (1 - p R ) + 2K?p R tanh 2 [2K? (1 - p R ) 



(32) 



(33) 



If K < K^, the disordered state M = will be stable, but if K > K^, a small perturbation 
will drive the system apart from the disordered state towards nonzero magnetization. Now 
we continue to find several interesting critical exponents. 

(1) x ~ \T — T C \~ J : Near the critical point, M eq — ► 0. If a weak field H is introduced, 
we will have 



fc 



ff (1 - Pfl) (Tfc+l + P/?T^T 2^ °j + 
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and 

dq k (t) 
dt 



-qk (t) + tann 

M 



K(l-p R ) (a k -i + a k+ x) + 2Kp R a + 



H 



k B T 



P({a};t). (34) 



Following the same way as that taken in the calculation of the critical point (in fact, we can 
just replace 2Kp R a by 2Kp R a + H/k R T), we get 
d 



dt 



M (t) = -M (t) + M (t) tanh [2K (1 - p R )] 



+ (2Kp R M (t) + JL) |l - 1 tanh 2 [2K (1 - p fl )]} . 



(35) 



From Eq . (|35|) one can easily find that, in a system in equilibrium near the critical point, 
K = K* + A, and 

dM 



X 



dH 



A" 1 . 



(36) 



Thus, 7 = 1. 

(2) M ~ |T - T/: When M -> 0, we take 

1 



AT 



J] (<7<7*_i<7 fc+ i) ~ M 3 . 



Then from Eq. (|3l| ) we can get 

« - (l - tanh [2# (1 - p R )\ - 2Kp R |l - ^ tanh 2 [2# (1 - p fl )]}) M (t) 
-i^tanh 2 [2K(l-p R )}M 3 (t). 

When K - Kf = A -»• 0, let dM (i) /dt = 0, and one will get M 2 ~ A by taking Taylor 
expansion. Thus, the critical exponent f3 = 1/2. 

(3) r ~ £ 2 ~ |T — Tcp 21 ": When studying the critical slowing down, we can assume M 
is very small and thus use Eq.(^). It yields 



M(t) = M(0)e- t/T , 



where 



r- 1 = 1 - tanh \2K (1 - p R )\ - 2Kp R \\ - ± tanh 2 \2K (1 - p R )}\ 
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If if -> iff, then r -> do. If If = iff - A, and A -> 0, then one will find r ~ A -1 . Thus 
= 1. It has been found in Monte Carlo simulations that v « 1/2, so z = 2. It is of the 

same value as that obtained on the regular lattice. 

ID Ising model on A-SWN show qualitatively the same behavior: Its critical point can 

be determined as 

tanh 2if c A + (N - 1) p A (l - ~ tanh 2 2K^j = 1, (37) 

and we have found the same critical exponents, 7, (3 and z. 

For a vertex on a A-SWN lattice, the number of the long range bonds is ha ~ (N — I) Pa, 
while for a vertex on a R-SWN lattice, n R ~ 2p R . The n- dependence of critical point if c 
of ID Ising model on A-SWN and R-SW, Eqs. fl33|) and fl37|), can be found in Fig. 1. 
One can clearly see the mentioned- above approximate relationship, oc — log 10 n A ~ 
-log M [(N - 1) p A ], or iff oc -\og w n R log 10 (2p R ). 

When ua and n R are small enough (in the small world region), from Eqs. fl33|) (|3?|) we 
can get 

n R ~ (l - tanh2iff) / Qiff tanh 2 2K? 



n A 



1 - tanh 2if c A ) / 



/v; 1 ( 1 -itanh 2 2if c A 



When and are approaching zero, if A ,B — > 00. For the same value of the critical point, 
(n R — ua) —> + . As shown in Fig. 1, for most of the region, the two curves are very close to 
each other. However, though the difference may be infinitesimal, the curve of n R is always 
above that of ha, as is distinct when they are relatively large. 



VI. THE INFLUENCE OF THE RANDOMNESS ON THE CRITICAL POINT 

On the behavior of the critical point, our results show interesting contrast between the 
Gaussian model and the Ising model. Here, we shall mention another interesting model 
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system, the Z)-dimensional mean-field (MF) Ising model, which might help us understand 
this problem. For a randomly selected spin in this model, each of its 2D nearest neighbors 
is replaced by an averaged one, and it is well known that the critical point is K c = 1/2D. 
A simple calculation will yield that on R-SWN will not change, while on A-SWN, since 
the long range bonds increase the contact of a spin with the system, K£ = 1/ {2D + tia)- 

Our result of the ID Ising model is in contrast to that of the MF Ising model and the 
Gaussian model, while for each one of them the critical temperature will be very close on 
A-SWN and R-SWN (the former will be higher), in the small- world region. This may be a 
result of the totally different role played by the long range bonds, (la) In the D-dimensional 
MF Ising model, the critical point is solely determined by the mean coordination number, 
which decides the coupling between an individual spin and the system. It is unchanged on 
R-SWN but will increase on A-SWN. Thus, the critical temperature will stay unaltered on 
R-SWN but will be slightly increased by the long range bonds on A-SWN, which typically 
take up only a small fraction. (16) As is shown by earlier studies the Gaussian model, 
though being a very different system, has a critical point which also only depends on the 
mean coordination number. On SWN, its critical point is very similar to the MF Ising model. 
It is believed to be mainly a mathematical result, and to understand this we shall review 
the calculations in Sec. |IV]. On R-SWN, for an individual spin, the long range interaction 
partly replaces the nearest-neighbor (n.n) coupling. However, in geography, the coordination 
number can be considered unchanged, since no bonds are created or eliminated (they are 
just redirected). As a result, on the right-hand side of the evolving equations of the spins, 
the lost part of the n.n coupling is exactly compensated by the MF term. The critical point 
is determined by taking average of the evolving equations of the spins, which only consist 
of linear terms. Thus, it is mathematically straightforward that the critical point will stay 
the same. On A-SWN, the consideration is similar, except that here the coordinate number 
will be slightly increased. (2) The cross-over observed in the one-dimensional Ising model 
is certainly governed by a different mechanism. For example, on R-SWN, there are two 
competing length scales: the correlation length £ ~ exp {2J/ksT), and the characteristic 
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length scale of the SWN, which can be taken as the typical distance between the ends of 
a shortcut |25|, £ ~ Vr^ D — Vr ■ When £ < (, the system basically behaves as a regular 
lattice, otherwise it will show MF behavior as an effect of the long range interactions. The 
transition occurs at £ « £, suggesting a critical point T c ~ |log 10 PH| _1 - Since the two 
structures, R-SWN and A-SWN, have the same length scales, the critical point can hardly 
be separated in the small-world region. However, interestingly we find that the mechanism 
in (1), i.e., long-range bonds — > coordination number (interaction energy) — > critical point, 
can also be observed, though being a minor factor. As is mentioned in Sec. [VQ , for the 
same expected number of long range bonds, the critical temperature on A-SWN will be 
higher than that obtained on R-SWN. We expect it to be a general phenomenon, though 
not yet reported by the numerical simulations. 



VII. SUMMARY AND DISCUSSIONS 

In this article, we study the critical dynamics, Glauber-type and Kawasaki- type, on two 
typical small- world networks, adding- type (A-SWN) and rewiring type (R-SWN). 

The logical sequence: As the introductory content we discuss the general approach of 
the critical dynamics, which is theoretically straightforward but may be mathematically 
too complex. We directly apply it to the kinetic Gaussian model governed by Glauber- 
type mechanism and obtain its evolution. We observe that, in the dynamic evolution of the 
individual spins, the influence of the system as a whole, which is the result of the presence of 
the long range bonds, takes the mean-field-like form as if it comes from an averaged spin. At 
the same time, earlier studies have revealed the mean-field (MF) static behavior of the Ising 
model and XY model. These both suggest us to present the following simplified method. All 
the SWN realizations are deemed as a single one, with both the effective Hamiltonian and 
the effective behavior averaged over all of them. It is tested in the same model and exactly 
leads to the same rigorous result. Then this method, which is believed to be theoretically 
reliable and mathematically feasible, is applied to two more difficult problems, the Gaussian 
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model governed by Kawasaki-type mechanism and the one- dimensional kinetic Ising model. 

The Gaussian model: (1) Whether it is built on A-SWN or R-SWN, the long range bonds 
introduce the mean-field-like global influence of the system to the dynamic evolution of the 
individual spins. On A-SWN such influence is additional while on R-SWN it partly replaces 
that of the neighboring spins. (2) On A-SWN, as more long range bonds are added, the 
critical temperature gets higher; but on R-SWN it does not differ at all from that of the 
regular lattice. This interesting discrepancy is explained in Sec. |VI[ (3) In both networks, 
the relaxation time, r = 1/ (1 — K/K c ), and thus zu = 1. This dynamic property has also 
been obtained on regular lattices and fractal lattices |T9| , pO| , [2T| , |22l . It is highly universal, 
independent of the geometric structure and the dynamic mechanism. (4) On SWN, the 
evolution of the Kawasaki-type model can be viewed as the combination of two mechanisms, 
the diffusion, and the automatic adjustment of the single spin to approach the average 
magnetization. The pure diffusion equation, -^q(r,t) = W 2 q (r,t), will be modified as, 
§- t q (r, t) = V'V 2 q (r, t) + C (M (t) -q(r,t)). By setting V and C to be zero we will get two 
competing characteristic temperatures, instead of the single definition of the critical point 
for the regular lattices (T> = 0). The temperature dependence of the dynamic evolution is 
discussed in Sec. |V B| . 



The Ising model analytically studied by the simplified method: The system shows very 
similar behavior on the two networks, A-SWN and R-SWN. Introducing very small per- 
turbation of the local magnetization to the disordered state, we obtain the critical point 
by judging the stability of the equilibrium. The inexistence of such a threshold p c for a 
finite temperature transition is confirmedB. From the dynamic equation we obtain the crit- 
ical exponents 7 and f3 in agreement with the numerical simulation and the already-proved 



x As mentioned above, for some physical considerations one may define a different J' for the long 
range bonds, which might be much smaller than J. With this method we can prove there is not 



such a threshold J' c either. 
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MF behavior of the system. The relaxation time is divergent near the critical point as 
t ~ \T — T c \~ l and thus zv — 1 (note the same relationship in the Gaussian model). In the 
ID Ising model, the SWN effect does not change the dynamic critical exponent z = 2. 

The influence of the randomness on the critical point: Our result of the ID Ising model 
is in contrast to that of the MF Ising model and the Gaussian model. For each one of them 
the critical temperature on A-SWN, T^, will be higher than that on R-SWN, T^, obtained 
for the same expected number of the long range bonds, though the difference may be hardly 
perceivable in the small-world region. A detailed analysis of the responsible mechanisms can 
be found in Sec. |VT[ 

Prospects: We hope the further studies of critical dynamics will continue to reveal in- 
teresting dynamic characteristics in the widely existing critical phenomena combined with 
SWN. The simplified method shall become a useful tool in this field. As the phase transition 
is of the MF nature, this method is certainly the best choice, and in this sense it is basically 
different from the custom MF approximation applied to other universality classes. Presently 
besides the numerical study, analytical treatment in the dynamic aspect is scarce compared 
to the study of the static properties |||| . In fact, at least in some cases, such a study may 
be feasible and fruitful indeed. Our work, especially that on the Ising model, also shows 
that the study of the dynamic aspect is often able to yield much information of the general 
properties, in a relatively convenient way. 
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APPENDIX A: 3D GAUSSIAN MODEL ON R-SWN 

(1). The two regular bonds on (connecting o~k-i and cr fe+1 ) have not been rewired, 
and the probability is (1 — p) 2 . This group can be further divided into many subgroups: 
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(1.0). There are no random bond on o k . The probability is (1 — p) 2 (l — -^-j- 



P 



N-2 



dq k K 

y- = ~Qk + y {Qk-i + Qk+i) ■ 

(1.1). There is only one random bond on a k . The probability is 
(1 - pf Clf^j^p (l - j^tp) N 3 , and 

dqk K 

-y = -qk + y {qk-i + qk+i + M) . 



(l.n). There is n random bonds 

2 / i \ n ( i \N—n—2 

on (Tjfc, and the probability is (1 — p) C^_ 2 [j^PI (1 — jf^iP) ■ This group can be 



further divided into C^_ 1 subgroups, each corresponding to n specific spins connected to a k 

2 / i \n / i \N—n—2 

via the random bonds, and with the same probability (1 — p) [j^ziPj (1 — j^rfP 
In the ith subgroup, 

dq k K( " \ 



Thus averaging them we get, 



-9fc + y (tffc-i + 9fc+i + nM ) 



l.(N — 2)). There is N — 2 random bonds on cr^, and the probability is 



:i-p) 2 GvM 



N-2 



-y = -9fc + y I 9fc-i + ^fc+i + E ft, I 

= -Qk + y + Qk+i + (N - 2) M) . 
Thus, the first part of the time derivative of single-spin, 



dt 



[1-pf 



K < 

-Qk + y (<?fc-i + qk+i) 



K 



N-2 



+1- E ^v- 



6 -o^ 2 ViV-l 



P 



1 - 



1 



TV - 1 



p 



N-n-2 



nM 



(Al) 
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(2). The bond connecting a k -\ an d a k is left unchanged but that connecting a k and 
a k+ i has been rewired. This bond can be rewired to each one of the N — 1 spins with equal 
probability (1 —p)p/N. In each case, we can analyze the situation in the way described 
above. For example, the bond is rewired to a y 

(2.j.0). There are no random bond on a k . The probability is 



N-2 



dq k K 



There is only one random bond on o k . 

N-3 



(2.J.1). 

sch^p (i - j^p) 

dq k 



dt 



K 

-Qk + y {Qk-i + Qj + M) 



The probability is 



(2.j.n). There is n random bonds on a k , and the probability is 

, N-n-2 



dq k K 

-y = ~qk + y (9fc-i + <7j + JiM) . 



(2.j.(N — 2)). There is iV — 2 random bonds on a k , and the probability is (j^p) 



N-2 



jQk = ~q k + j (q k -i + qj + (N - 2) M) 



Thus the second part of the time derivative of single-spin, 



K 



N-2 



n 

u n=0 



^ v N-n-2 

-p | nM 



(A2) 



M-U V iV-r 

(3). The bond connecting and o"fe+i is left unchanged but that connecting a k -i and 
a k has been rewired (we omit the very small probability that this bond may be "rewired" 
to a k ). Now there is only one regular bond on a k . (Pay attention that this case is different 
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from (2).) Based on similar consideration, we have the third part of the time derivative of 
single-spin obeys 



dq k 
dt 



[1 -p)p 



K 

-Qk + 



K N - 2 

U n=0 



P 



N-l 



N-l 



P 



N-n-2 



nM 



(A3) 



(4). Both the bond connecting ak&k+i and <Jk-i&k have been rewired. Based on similar 
consideration, we have the fourth part of the time derivative of single-spin 



' dg k 
dt 



V 



K K N ~ l 



n=0 



V 



N-l 



1 - 



N- 1 



V 



N-n-2 



nM 



(A4) 



Applying Eq. fllCf) we get 



dqk 
dt 



E 



dq k \ 
dt I . 



2K »r „ sK , 

-Qk + P~rM + (l-p) — (g fc _! + q k+l ) 
b b 



(A5) 
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Caption of figures 

Fig.l. The n-dependence of critical point K c of ID Ising model on A-SWN and R-SWN. 
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